Magnetoelastic coupling in triangular lattice antiferromagnet CuCrS2 
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CuCrS2 is a triangular lattice Heisenberg antiferromagnet with a rhombohedral crystal structure. 
We report on neutron and synchrotron powder diffraction results which reveal a monoclinic lattice 
distortion at the magnetic transition and verify a magnetoelastic coupling. CuCrS2 is therefore an 
interesting material to study the influence of magnetism on the relief of geometrical frustration. 
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INTRODUCTION 



Heisenberg antiferromagnets on a triangular lattice 
are subject to ongoing interest due to a strong cor- 
relation between lattice geometry, electronic and mag- 
netic properties. Geometrical frustration plays a key 
role in compounds with triangular arrangement of mag- 
netic moments^ and leads to particular characteristics 
including incommensurate spin structures and multifer- 
roic properties^ Ternary triangular lattice dioxides and 
dichalcogenides have in common to crystallize in a lay- 
ered structure with strong crystalline anisotropy perpen- 
dicular to the layers (mainly space group R3m and R3m) . 
In chromium based dioxides the most commonly estab- 
lished magnetic structure is a quasi two-dimensional 120° 
spin structure, typical for Heisenberg exchange on a tri- 
angular lattice, with a weak interlayer coupling e.g. in 
AgCrOA PdCrOA NaCr0 2 5 and CuCrOA The latter 
two also show a magnetoelectric coupling where the spin 
structure induces ferroelectricity 7 ^. Ternary chromium 
dichalcogenides have been less intensively studied, and 
show a variety of different magnetic structures ranging 
from a 120° spin structure in LiCrS^ over a commen- 
surate magnetic structure in KCrS^ to a helix with an 
in-plane spin orientation in NaCrS2^°- 
The compound under investigation, CuCrS2, has drawn 
attention as an ionic conductor 10 and is assumed 
to exhibit a spin glass state under the substitution 
of chromium ions by vanadium i 11 ' 12 A microscopic 
temperature-dependent investigation of magnetic and 
crystallographic properties is, however, lacking. CuCrS2 
crystallizes in space group R3m at room temperature^ 
The basic atomic structure consists of covalent-ionic 
bound S-Cr-S layers separated by a van der Waals gap 
causing strong crystalline anisotropy perpendicular to 
the layers. Isolated CrS2 is metastable and occurs 
only with an electron donor, therefore monovalent Cu + 
cations are intercalated between (CrS2 )-sandwiches . 10 i 14 
The magnetic Cr 3+ ions with spin S = 3/2 form trian- 
gular layers which are shifted by (|, |, |) and (§13,3), 
as shown in the inset of Fig. [TJ The nearest-neighbor 
intralayer Cr-Cr distance of 3.48 A is small compared 



to the interlayer Cr-Cr distance of 6.55 A. The mag- 
netic exchange within the layers over Cr-S-Cr pathways 
is expected to be much stronger than the interlayer ex- 
change via Cr-S-Cu-S-Cr bonds. Hence, one would ex- 
pect to find the crystalline anisotropy to be reflected in 
the magnetic system, as described in the dioxides. Never- 
theless, early neutron powder diffraction experiments on 
CuCrS2 revealed a three-dimensional magnetic ordering 
below TV = 40 K into a complex helical structure with 
an incommensurate magnetic propagation vector.— This 
phenomenon cannot be explained without interlayer ex- 
change interactions which must be of the same order of 
magnitude as the intralayer exchange. An estimation of 
exchange constants is given in the discussion. The partic- 
ularity of CuCrS2 is a strong lattice effect exactly at Tjy. 
We assume that the lattice distortion provides a chan- 
nel for the system to escape the 2D triangular arrange- 
ment, usually leading to magnetic frustration. The conse- 
quences on the lattice and magnetic structure throughout 
the phase transition are subject of this paper. A similar 
lattice effect occurs in triangular lattice CuFe02 , where a 
crystal symmetry lowering is found to lift the degeneracy 
of the frustrated spin system 1 ^ together with multiferroic 
properties induced by a magnetic field or non magnetic 



impurities 
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Strong spin-lattice effects are also found 



in 3D frustrated systems, as shown for the well-known py- 
rochlorc spinel ZnCr204. 19 The compound was found to 
undergo a Spin-Peierls-like phase transition which com- 
prises an energy lowering through a lattice distortion and 
the opening of an energy gap. The origin of this gap 
has been explained by clustering of the Cr moments into 
weakly interacting antiferromagnetic loops^S Although 
the formation of spin cluster has been clarified, the mag- 
netic structure is not fully understood^ We claim to see 
a similar effect on a triangular lattice and report in this 
paper on the detailed magnetic structure. 



II. EXPERIMENTAL DETAILS 

The polycrystalline sample of CuCrS2 was synthesized 
by solid state reaction of the pure elements. A de- 
tailed description is given elsewhere^. Single crystals 
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FIG. 1: (Color online) Magnetic susceptibility of CuCrS2 
measured at H = 10 kOe as a function of temperature. An 
antiferromagnetic anomaly at Tjv = 37.5 K is clearly visible. 
The red solid line shows the Curie- Weiss fit to the paramag- 
netic part of the curve. Inset: Structural unit cell of CuCrS2 
at room temperature with space group R3m. 



were grown by chemical vapor transport in an evacuated 
quartz ampoule with p < 10 mbar. A temperature 
gradient was applied over 30 cm ampoule length between 
900°C and 950°C during three weeks. The furnace was 
then cooled at a cooling rate of 1 K/min to room tem- 
perature. The resulting single crystals are thin platelets 
which possess a black shiny surface with the [001] direc- 
tion perpendicular to it. The size of the surface ranges 
from 10 to 100 mm 2 with a typical thickness of 0.2 mm. 
The stoichiometry of the samples was checked by energy 
dispersive X-ray spectroscopy and verified in case of the 
polycrystalline powder. The single crystals contain an 
impurity phase of not more than 10% CUCT2S4. System- 
atic variations of growth conditions in order to avoid the 
impurity phase, produced crystals too thin (< 10 /im) 
and unstable for use in neutron scattering experiments. 
Neutron powder diffraction measurements were carried 
out on the cold diffractometer DMC at the neutron spal- 
lation source SINQ, Paul Scherrer Institut (PSI) and 
the thermal diffractometer DlA at the Institut Laue- 
Langevin (ILL) with A = 2.46 A and A = 1.91 A , re- 
spectively, in a standard orange cryostat. Single crys- 
tal diffraction experiments were performed on the ther- 
mal instrument TriCS^ (SINQ) with A = 1.18 A in a 
four-circle Eulerian cradle and on the three-axis spec- 
trometer IN3 (ILL) with A = 2.36 A. High-resolution 
synchrotron radiation powder diffraction patterns were 
collected at the powder diffraction station of the Swiss 
Light Source Materials Science (SLS-MS) beamline at an 
incident photon energy of A = 0.621288(1) A. The sample 
was mounted in a Janis cryostat in 0.2 mm Lindemann 
capillaries which spun at approximately 10 Hz during 
the 28 scan to avoid preferred orientation. The collected 
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FIG. 2: (Color online) Powder diffraction data taken on DlA 
at 10 K (refined) and 44 K (raw data). The refinement of 
the 10 K data with a nuclear and magnetic phase is indicated 
by a black solid line. The difference between refinement and 
observed data at 10 K is shown in grey, (a) depicts the low 29 
region where magnetic reflections are visible which disappear 
at 44 K and (b) shows the high 29 region where a splitting of 
the nuclear Bragg peaks below the transition is observable. 

neutron and X-ray diffraction data were refined using the 
program Fullprof 24 based on Rietveld refinement. 
Magnetic susceptibility measurements were carried out 
on a physical properties measurement system over a tem- 
perature range of 4 — 300 K with an external field of 
10 kOe and on a SQUID magnetometer for temperatures 
up to 400 K. 

III. RESULTS 

Figure [T] shows the temperature dependence of the 
magnetic susceptibility of the powder sample with an an- 
tiferromagnetic transition at T/y = 37.5 K. From the lin- 
ear part of the curve the paramagnetic Curie- Weiss tem- 
perature could be determined to &cw — —118.5(1) K 
and the effective magnetic moment amounts to [i eX p — 
3.75(1) [Xb close to the spin-only value /is 3/2 = 3.87 \xb- 
The magnetisation curve of the single crystal, mea- 
sured up to 400 K (not shown), indicates a superpo- 
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FIG. 3: (Color online) (a) Correlation of the rhombohedral 
unit cell in hexagonal description (R3m) and the monoclinic 
unit cell (Cm), (b) Influence of the monoclinic angle (3 on the 
layers. The angle is overdrawn to show the effect of distortion; 
actually /3 only deviates slightly from 90°. 



sition of the antiferromagnetic transition of CuCrS2 at 
TV = 37.5 K and a ferromagnetic transition of CuCr2S4 
at T c = 377 

A good measure for the degree of frustration is the ra- 
tio \®cw\/Tn which is in case of CuCrS2 only 3.16 and 
therefore rather small to infer a strongly frustrated sys- 
tem with values usually larger than 10i2& 
Neutron powder diffraction data taken on D1A and DMC 
at different temperatures show strong magnetic Bragg 
peaks emerging below the magnetic transition, i.e. Tjy = 
37.5 K. This observation must be related to a magnetic 
ordering of the system. At higher diffraction angles 29, 
accessible on D1A, a splitting of nuclear Bragg peaks is 
also observable, which indicates a symmetry lowering of 
the crystal structure. The pattern taken above and below 
the transition at 44 K and 10 K, respectively, are shown 
in Fig. O The pictures (a) and (b) show the refinement 
of the low temperature data at 10 K with two phases, nu- 
clear and magnetic, superimposed by the raw data pat- 
tern at 44 K. Above TV the Cr 3+ moments show no long 
range order anymore, but localized Cr 3+ moments still 
cause paramagnetic scattering. Therefore, the difference 
in background intensity between spectra taken at 10 K 
and 44 K at low angles 29 (Fig. [DJa)) can be explained 
by an angle dependence of the paramagnetic scattering 
due to the magnetic form factor. The effect is negligible 
at higher angles 29 (Fig. Mb))- 

The structural transition indicated by a splitting of the 
nuclear Bragg peaks in the low temperature phase could 
be identified as a symmetry lowering from rhombohe- 
dral space group R3m (160) to monoclinic space group 
Cm (8). The relation between both unit cells is given in 
Fig. [31 as well as a visualization of the monoclinic angle (3. 
The introduction of a monoclinic angle (3 causes a slight 
shear movement of the layers in the afr-plane. Lattice pa- 
rameters and crystal symmetry at 300 K and 10 K taken 
at the SLS-MS beamline are listed in Table |TJ Refine- 
ments with space group Cm describe the splitting of the 
Bragg peaks in the nuclear phase very well (see Fig.^b)). 
Above the phase transition all three lattice parameters, 
depicted in Figs. BJa) to[2|c), decrease with tempera- 




FIG. 4: (Color online) Temperature evolution of the lattice 
parameters a, b and c corresponding to (a) to (c), extracted 
from synchrotron data, in monoclinic description Cm. For 
better comparability the refinements for T > 37 K were 
performed in the monoclinic description with j3 = 90° and 
a — -J?,b. (d) shows the refined weight percentage of both 
the high temperature R3m phase and the low temperature 
Cm phase through the transition. The error bars are smaller 
than the symbols. 

ture. Upon further cooling a and c show a discontinuous 
jump towards lower values, whereas b increases quickly. 
The complementarity of the in-plane lattice constants a 
and b will later be explained by a distortion of the Cr 
hexagon. 

TABLE I: Lattice parameters of CuCrS2, refined with space 
group R3m at 300 K and Cm at 10 K. Data were collected 
at the SLS-MS beamline with wavelength A = 0.621288(1) A. 



Temperature 


300 K 


10 K 


Space group 


R3m 


Cm 


Nr. 


160 


8 




rhombohedral 


monoclinic 


a [A] 


3.484543(4) 


5.99923(1) 


b[A] 


3.484543(4) 


3.491502(5) 


c[A] 


18.71878(2) 


18.64794(3) 


a [°] 


90 


90 




90 


89.93727(8) 


7[°] 


120 


90 


RBragg[%] 


4.13 


3.83 



The refinement of magnetic intensities was done with a 
magnetic propagation vector documented earlier— which 
was redefined on the monoclinic lattice and fits well to 
the magnetic Bragg peaks (see Fig. Ufa)). The propaga- 
tion vector k= (-0.493,-0.087,1.25) in relative length 
units (r.l.u.), describes a three-dimensional helical ar- 
rangement of magnetic moments originating from Cr 3+ 
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FIG. 5: (Color online) Intensity of the strongest magnetic 
Bragg peak at 29 = 13.96° (squares) and the monoclinic an- 
gle j3 (circles) as a function of temperature, extracted from 
neutron powder diffraction data taken on D1A. While cooling 
the sample, the ordering of magnetic moments and the lattice 
distortion set in simultaneously. 



ions. The best fitting results were obtained with no 
magnetic moment on the Cu atoms, from which we con- 
clude that Cu is in oxidation state Cu + , corresponding 
to the stoichiometry of the sample. Neutron single crys- 
tal diffraction experiments were performed on IN3 (ILL) 
and TriCS (SINQ) to confirm the rather unusual mag- 
netic propagation vector and to avoid artifacts from the 
powder refinement. Magnetic Bragg peaks were found at 
several positions in 3D reciprocal space, corresponding to 
the propagation vector obtained from powder measure- 
ments and verified the incommensurate magnetic struc- 
ture. The intensity of all magnetic peaks, including the 
strongest (003)-fc, vanishes at Tjv = 37.5 K but shows 
different temperature behavior than in the powder, i.e. a 
varying critical exponent extracted from M s ~ y/7 close 
to the phase transition. The values are 0.08(4) for the 
powder compared to 0.5(1) for the single crystal. This 
is possibly due to an internal magnetic field in the single 
crystal caused by the ferromagnetic impurity CuCr2S4 
which slows down the transition. The magnetic inten- 
sity is saturated for T < 30 K. Since TV and k coin- 
cide for powder and single crystal in both susceptibility 
and diffraction data, we conclude that the presence of 
CuCr 2 S4 does not disturb the measurement of the mag- 
netic structure in CuCrS2 at low temperatures. 
A detailed analysis of the DlA powder data revealed 
an interesting coupling between magnetic ordering and 
structural distortion. The temperature dependencies of 
the magnetic intensity in the strongest Bragg peak at 
29 = 13.96° and the monoclinic angle f3 show reverse 
behavior, depicted in Fig. [5j As the sample is cooled 
below Tn the magnetic ordering sets in and a deforma- 
tion of the lattice develops simultaneously, expressed by 
a monoclinic angle f3 deviating from 90°. Although the 
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FIG. 6: (Color online) Temperature dependence of in-plane 
Cr-Cr bonds measured by synchrotron powder diffraction. In- 
set: Different kinds of chromium bonds and their influence on 
the distortion. The deformation is overdrawn for better illus- 
tration. 

deviation of (3 from a higher symmetric value seems to 
be tiny it is crucial to solve the low temperature crystal 
structure. Figure [5] shows the mutual influence of lattice 
distortion and magnetic ordering. 

Additional to the interlayer shear movement, another 
type of intralayer structural distortion was identified 
from synchrotron powder measurements at the SLS-MS 
beamline. At high temperatures the in-plane chromium 
lattice describes a symmetric hexagon. Upon cooling 
the symmetric chromium environment distorts, such that 
four of the six equal Cr-Cr bonds contract and two elon- 
gate. This causes a flattening of the hexagon, as schemat- 
ically shown in the inset of Fig. [6j and leads to nearest- 
neighbor (di) and next nearest-neighbor (cfe) Cr-bonds. 
The Cr-Cr distances d\ and g?2 show the same tempera- 
ture evolution as the interlayer distortion (see Fig. [6]). In 
addition, di coincides with the lattice parameter b, which 
explains its unconventional temperature dependence be- 
low Tn. 

These two types of distortions, the interlayer shear 
movement by an angle and the intralayer hexagonal 
symmetry breaking, seem to play an important role for 
the magnetic ordering. The system is able to select a 
magnetic ground state due to the relief of geometrical 
frustration. 



IV. DISCUSSION 

The analysis of the diffraction pattern at different tem- 
peratures shows an abrupt change of Cr-Cu (lattice pa- 
rameter c) and Cr-Cr distances (da) below the phase tran- 
sition. Usually a contraction causes strain in the lattice 
which is proportional to the square of the atomic dis- 
placement u and increases elastic energy ~ u 2 . A phase 
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transition has to be energetically favorable, which means 
that the gain in elastic energy must be compensated. Us- 
ing the example of the phase transition in solid oxygen, 
Rastelli et al£L state that a first-order transition is fea- 
sible when the ground-state energy of the rhombohedral 
phase equals the ground-state energy of the monoclinic 
phase. Synchrotron powder diffraction data confirm a 
first-order transition in CuCrS2 identified from a resid- 
ual rhombohedral Rim phase even at the lowest tem- 
perature, depicted in Fig. Q] (d). Mixed-phase regimes 
are characteristic for first-order phase transitions^. This 
may explain the deviation of the critical exponent 0.08(4) 
from common models (Ising, XY or Heisenberg), since it 
is merely defined for second-order phase transitions. A 
reduction of the total free energy can only be achieved 
by a decrease in magnetic energy. For discussion we as- 
sume three different magnetic exchange paths in CuCrS2 
comprising in-plane nearest-neighbors J\ (along d\, see 
Fig. and next nearest-neighbors J 2 (along (I2) and 
inter-plane exchange J3. 

The intralayer exchange integrals J\ and J2 can have two 
origins, a direct antiferromagnetic interaction of neigh- 
boring Cr 3+ ions or an indirect superexchange interaction 
via S 2_ which would be ferromagnetic due to the Cr-S- 
Cr angle of 90°. The susceptibility curve shows antifer- 
romagnetic behavior, and we can therefore neglect the 
ferromagnetic superexchange path, fnterlayer exchange 
is presumably mediated by Cu + ions through the path 
Cr-S-Cu-S-Cr. It needs to be taken into account to ex- 
plain the three-dimensional magnetic structure. The con- 
traction of Cr-Cu bonds through the transition seems to 
provide evidence for the important role of the interlayer 
Cr-Cu-Cr exchange path. The ratio between the in-plane 
and intcr-planc exchange constants can be calculated for 
a rhombohedral antiferromagnet with a helical propaga- 
tion vector £1 For CU&S2 we find a ratio J3/J1 = 2, 
which would imply a very strong inter-plane coupling, 
but this value needs to be confirmed by inelastic neutron 
scattering. 

The in-plane distortion of the high temperature symmet- 
ric Cr 3+ -hexagon consists of a deformation where four 
nearest-neighbors move towards the central Cr 3+ ion and 



two nearest-neighbors move away. The transition takes 
place in order to maximize magnetic energy, as was found 
in solid oxygen^. The two different kinds of intralayer 
Cr-Cr distances d\ and c?2 (Fig. may be evidence for 
the formation of magnetic clusters in CuCrS2, which is 
supported by preliminary neutron scattering results. It 
has been found that localized magnetic excitations occur 
additionally to standard spin waves branches in the spec- 
tra. To our knowledge the formation of magnetic clus- 
ters in triangular antiferromagnets with spin-lattice cou- 
pling has only been theoretically described™ so far. This 
phenomenon makes CuCrS2 an interesting candidate for 
studying new magnetic effects on frustrated lattices. 
V. CONCLUSION 

From the quasi two-dimensional layered structure of 
CuCrS2 one would expect a low dimensional magnetic 
ordering due to the van der Waals gap between magnetic 
Cr 3+ layers. Nevertheless, our measurements show un- 
ambiguously that magnetic moments in CuCrS2 order in 
a fully three-dimensional manner with a helical magnetic 
propagation vector fc= (—0.493,-0.087,1.25). Moreover, 
the system undergoes a crystallographic transition from 
rhombohedral Rim to monoclinic Cm at Tjv = 37.5 K, 
which coincides with susceptibility measurements. Since 
the ordering of magnetic moments and the lattice distor- 
tion occur simultaneously, CuCrS2 has been identified as 
magnetoelastic material. The system is therefore an in- 
teresting candidate for studying spin-lattice effects on a 
triangular lattice. 
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